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Yucesoy, Katzgraber, and Machta Reply: Billoire et al. lHJ 
criticize the conclusions of our Letter t2J. They argue that 
considering the Edwards-Anderson (EA) and Sherrington- 
Kirkpatrick (SK) models at the same temperature T is inap- 
propriate and propose an interpretation based on replica sym- 
metry breaking (RSB). In our Letter we compare the SK and 
EA models at the same low reduced temperature T « QATc- 
Billoire et al. compare them at different T such that P{q = 0) 
is nearly the same. They also consider the quantity A{qQ, k) 
12], which measures the probability with respect to the dis- 
tribution of couplings ,/ that Pj{q) exceeds k in the range 
l^l < go- In the low-T phase A ^ if a two-state picture 
holds, while A ^ 1 if RSB holds. Considering the same T 
for both models was not essential to our argument; however, 
we think it is important to study both models at the lowest 
temperature possible to understand the low-T phase. 

For the EA model, we simulated systems up to size L = 12 
at r = 0.423, whereas Billoire et al. simulated sizes up to 
L = 32 but at T = 0.703. We find A leveHng off as a function 
of L at low T (see Fig. 5, Ref. t2J); they find it increasing as a 
function of L at higher T (Fig. 1 inset, Ref. [T]). It is not clear 
which trade-off in L vs T better reflects the infinite-volume 
behavior. However, P{q) for L = 12 and T — 0.42 is closer 
to a (5 function at ^ea, which is the infinite-volume behavior: 
P((7ea) divided by the width at half maximum of the qea 
peak equals 29.1 for L = 12 at T = 0.423, and 18.4 for 
L = 32 and T = 0.703 ID. We also note that the increase in 
A seen in Fig. 1 (inset) of Ref |T| is most pronounced for L = 
32. However, this point appears to be anomalous and P{q) 
from the same simulations [3 1 shows a similar anomaly, which 
may reflect large statistical errors or incomplete equilibration. 
Finally, Ref. UJ studies bimodal disorder, which converges 
slower m than Gaussian. 

The theory in the Comment attributes the plateau in A for 
our EA data to a combination of a small value of /(go) — 
J\ |< P{<i)dq and the slow growth in L of F((7ea)- It pre- 
dicts that A for the EA model will grow to unity extremely 
slowly in L. Our Fig. 6 shows that even after factoring out 
the slower growth in P{q-EA) for the EA model compared to 
the SK model, we still find a qualitative difference between 
the two. The proposed RSB scaling theory 1 1 1 asserts that 
A{qo,K) - [P(gEA)/K]^^*' - 1- This can be simplified 
when I{qo) is small. Noting that /(go) ~ QqP{0) one ob- 
tains A(go,K) « go/^(0)log[F(gEA)/K]- The predicted lin- 
ear dependence of A on go is consistent with our data and is 
neither surprising nor a strong test of the theory. The fact that 


data from different sizes lie on similar curves agrees with the 
plateau in our data but does not demonstrate that A is actually 
growing slowly with L for fixed go and k. To test the validity 
of the proposed theory 1 1 1, we compared our data for the EA 
model at several T holding P(gEA)/K and /(go) fixed. The 
proposed theory predicts that if these variables are fixed, A 
should remain constant. Figure [T] shows A(go,K) vs T for 
L = 10, 12. For each T and L, both go and k are adjusted so 
that /(go) ~ 0.067 and P{qEA)/K — 3 (go ranges from 0.16 
to 0.56 and k from 0.5 to 2.6 as T decreases from 0.7 to 0.2). 
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FIG. 1: (Color online) A(go,«:)vsr for L = 10 and 12 (EA model). 
Qo and K are chosen such that /(go) ~ 0.067 and P(gEA)/K — 3. 


Figure [TJshows that A is not constant as predicted by Ref. 1 1 J, 
and therefore this theory does not explain our EA data. 

In conclusion, we stand by our assertion that the low- 
temperature behavior of the EA model does not appear to be 
mean-field-like. 

B. Yucesoy,^ H. G. Katzgraber,^ and J. Machta^ 

^Physics Department, University of Massachusetts, Amherst, 
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